Topology and Energetics of Metal-Encapsulating Si Fullerene-Like 

Cage Clusters 

o, 

O [ Takehide Miyazaki", Hidefumi Hiura^ and Toshihiko Kanayama" 

W). ""Advanced Semiconductor Research Center, National Institute of Advanced Industrial Science and 

■^ , Technology, AIST Tsukuba Central 4, 

a^. 

1-1-1 Higashi, Tsukuba 305-8562, Japan 

^ , Fundamental Research Laboratories, NEC Corporation, 

34 Miyukigaoka, Tsukuba 305-8501, Japan 

(February 1, 2008) 
-)— > ■ 

s 

jIh ■ Abstract 

O 

o 



a 



On the basis of a topological discussion as well as an ab initio calculation, 
we show that it is possible to construct a fullerene-like Si cage by doping of 



(N 
> 

l/^ , a transition metal atom. The cage is a simple 3-polytope which maximizes 

o 



the number of its inner diagonals close to the metal atom. Our topological 
argument also reveals how closely the structure of the fullerene-like Si cages 
studied is related to that of fullerenes themselves. 
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Synthesis of fullerene-cage-like clusters composed of elements other than carbon (C) is a 
subject of great interest ^. Especially, it is intriguing to study whether silicon (Si) analogues 
of the fullerenes can exist in energetically favorable forms. Production of Si cage clusters is 
also important from a technological point of view. They may be used as building blocks for 
fabrication of various structures in electronic devices in ten nanometer scales. 

A main question which we would like to answer in this study is how one can obtain a 
smooth fullerene-like cage composed of Si atoms. Due to their sp^ nature, the Si atoms tend 
to bind themselves against generation of fullerene cages. As a matter of fact, Si„ clusters (n 
up to ~50) usually favor compact forms which are completely different from fullerene cages 
0J^. The structure of a Sieo cluster in a fullerene cage is highly distorted 0]. 

As is well known, the structure of a typical fullerene cage (C„) is composed of only 5- 
and 6-membered rings. The number of the 5-membered rings is twelve irrespective of n. 
Each atom is connected with three neighbors. Stability of fullerenes may be attributed to 
two factors: a bonding among sp?, hybrids, where p\\ stands for the intra-cage component of 
the C p orbital, and vr conjugation among p-orbital components normal to the cage surface 

iP±)- 

A fundamental difficulty to maintain a smooth Si fullerene cage stems from the fact that 
the TT conjugation among the p± orbitals does not occur. Thus, the caging mechanism is 
only the a bonding of the spy hybrids. The distortion of the fullerene cage of Sieo is due to 
admixture of substantial p± components with spf,. 

A promising solution to the augmentation of Si cage structures is to put one or more 
additional atoms inside of the cages. If electron orbitals of these "doped" atoms have 
a substantial overlap with the p± orbitals, then the additional cohesion forces would be 
supplied to the cage. This idea has originally sprouted from the construction of Si„ clusters 
with n larger than ~20 to account for the exceptional chemical inertness of Siag and Si45 
01 and also for "prolate-oblate" structural change in Si„ around n^27 . A common aspect 
of the Si„ configurations studied in Refs. is that they are configured to mimic as much of 
the bulk Si structure as possible, while the outermost cages yet resemble the counterparts 
of fullerenes. 

As explained above, there are at least two mechanisms to stabilize a doped Si fuUerene- 
cage cluster: (a) a bonding within the sp?, hybrid network of the cage and (b) bonding 
among the pj_ orbitals of the cage and the orbitals of doped atoms. A key issue to be solved 
for production of stable fullerene-like Si cage clusters is how these two factors should be 
tuned in order to maximize the total binding among the constituent atoms in the clusters. 

The goal of this study is to give an explicit answer to this question from a theoretical 
point of view. For this purpose, we target metal-doped Si cage clusters, M@Si„, because 
we have recently experimentally suggested that stable metal-encapsulating Si cage clusters 
(M@Si„, M=Hf, Ta, W, Re, Ir, etc. and n<~14) might be grown g. We also note the 
outstanding property of M@Si„. When doped with Na and Ba, the M@Si„ clusters in 
certain phases of the solid state exhibit superconductivity [Q. Although some theoretical 
T0| as well as experimental |[ll| studies have been performed for metal-doped Si clusters. 



roles of the cage topologies played in stabilization of fullerene-like Si cages have been little 
argued. 

First, we discuss what topology should be suitable for a cage of M@Si„ (n<20). A 
very important clue to the solution has been recently given by mathematicians, Bremner 



and Klee |]T2[. They study inner diagonals (see below or Ref. ||T2[ for definition of an inner 



diagonal) of convex polytopes. The quantity of our own interest here is the number of inner 
diagonals denoted by 63. We will propose that a simple 3-polytope [jl3[, in which the number 



of inner diagonals close to the dopant (^3 ) is maximized, should be a candidate of a stable 
fuUerene-like Si cage of M@Sin. 

Next, we perform a first-principles energetics of the clusters in the case of M =W (W@Sin, 
8<n<16). The clusters in fuUerene-like cages predicted by the above topological picture 
are indeed obtained for n=12 [Fig.l(d)] and n=14 [Fig.l(e)]. For ?7,=8 and 10, however, the 
caging with corresponding simple 3-polytopes [Fig. 1(a) and Fig. 1(b)] may not occur, because 
their total energies are significantly higher than those of the respective isolated atoms. For 
n=16, a single cage is unstable (Fig. 2). The cage is relaxed into a skeleton similar to either 
of n= 12 [Fig. 1(d)] or n=14 [Fig.l(e)]. Thus we conclude that it is possible to construct 
energetically favorable fuUerene-like Si cage clusters doped with a transition metal atom, 
whose cage sizes depend on the dopant. 

Now we describe our results. First a topological discussion of Si cages is given. We 
denote the number of p-membered rings of a Si cage as Np. Throughout this study, we 
limit the range of p to 4<p<6. We also use a vector notation N={N4,N5,Nq) for brevity. 
For example, simple 3-polytopes with n=8 and 10 vertices are a cube [iV=(6,0,0), Fig.l(a) 
|T^] and a pentagonal prism [iV=(5,2,0), Fig. 1(b)], respectively. For each n>12, simple 



3-polytopes with different A^'s exist as shown in Fig.l and Table I |T5 



An inner diagonal (or a 3-diagonal) of a convex 3-polytope P is defined as a segment 



that joins two vertices of P and that lies, except for its ends, in P's relative interior [|12 
There are two other kinds of diagonals of P, 1- and 2-diagonals. The former are the edges. 
As for the latter, there are two, five and nine 2-diagonals in a 4-, 5- and 6-membered ring, 
respectively. If we denote the number of i-diagonals by 6i, then 61+52+6^= 2 ' where 
n is the number of vertices. For example, 6^=4 for A^=(6,0,0) and 6^=10 for A^=(5,2,0). 
In this study, we further define 6^ as 6^ =J2i>jd{dcut — dij), where i and j are indices 
of Si atoms, dent is the cut-off radius and dij is the distance between a metal atom and a 
diagonal joining the i-th and j-th Si atoms. 6{x) is a step function; 6{x)=l for x>0 and 
6{x)=0 otherwise. Obviously, 6^-'^ approaches "'^"~ -* for dcut^>-oo. In order to count the 
number of only inner diagonals, dcut niust be small. In practice, we vary dcut from I.OA 
down to 0.25A to check the dcut-dependence on 6^ . The overlap between the p± orbitals 
of a Si cage and the counterparts of the dopant metal atom may be roughly proportional 
to the number of the p± orbitals pointing toward the metal atom, which should also be 
approximately proportional to 6^^^ at a small dcut by its definition. Thus analyzing 6^-^^ 
may measure the strength of the bonding between the doped atom and the Si cage. Our 
first-principles calculation shows that 6^ at 5cut=0.25A is maximized for the lowest-energy 
W@Si„ cluster for each of n=12 and 14. 

Bremner and Klee |]12[ show that, for a given n, P with a maximum ^3 is simplicial 
13| . A simplicial 3-polytope seems the best for a cage of M@Sin for each n, because the 



energy gain due to the factor (b) should be fully enjoyed. This is not true, however, since 
the mechanism (a) is instead sacrificed so that a total energy gain is not maximized. Our 
first-principles calculation shows that a W@Sii2 cluster with an icosahedron cage {i.e. a 
simplicial 3-polytope with 12 vertices and ^3=36) is energetically very unfavorable [5.76 
eV higher in energy than the (6,0,2)-cage cluster shown in Fig. 1(d)], despite that all p± 



orbitals of the cage point toward the central W atom. The averaged Si-Si distance (^2) of 
the icosahedron cage (2.72A) is much longer than those of both (4,4,0)- {d2=2A2A) and 
(6,0,2)-cages {d2=2.39A). In contrast, the averaged Si-W distances (di) are similar, ~2.6A. 
This evidences that too much preference of the factor (b) is taken over (a) in the icosahedron 
cage. Therefore a simple 3-polytope is a good candidate for a Si^ (n<20) cage, in which 
both (a) and (b) factors may work equally well to enhance the stability of the M @Si„ cluster 
as a whole. 

Turning to the topology of a carbon fuUerene cage, its peculiarity can be characterized 
in light of 63. Let Q be a simple 3-polytope with n vertices. For each of n>24, Q with F 
(= n/2 + 2) facets composed of only twelve pentagons and F— 12 hexagons maximizes ^3 
which is n{n — 13)/2 -|- 30 |0. For n=20, Q with the largest 63 (=100) is a dodecahedron 
which has only twelve pentagons. In the case of n=22, however, Q with A^5=12 cannot be 
realized [|12|. The maximum possible value of 63 is 128 when A^=(l,10,2). Thus a typical 



fuUerene C„ cage (n>20) can be identified to be a simple 3-polytope with a maximum value 
of ^3 except n=22. The Si„ cages we consider for 8<ri<16 include Q's with maximum 63 s. 
Although 63 (=18) of the (6,0,2)-cage [Fig.l(d)] is not a maximum, it is the second largest 
value. This is why we can say that our Si cages are fullerene-like. 

In order to substantiate the above topological argument, we perform a first-principles 
energetics of the W@Si„ clusters. Here we adopt a recently proposed single-parent evo- 
lution algorithm (SPEA) [|l7l to update the atomic coordinates in an unbiased way. For 
SPEA simulation, we calculate the total energies of clusters with a linear combination of 
atomic orbitals (LCAO) using the GAUSSIAN98 package [l^. Then we use the calculated 
coordinates of the clusters as the inputs of plane- wave (PW [|I^) total energy calculations 
20| where quenched molecular dynamics is performed for final convergence. Electronic 



structures of clusters are calculated with density-functional theory [^ within generalized 
gradient approximation to the exchange-correlation energy functionals (Becke'88 [^ and 
Perdew-Wang'91 [g§ for LCAO and Perdew-Wang-Ernzerhof'96 Q for PW) |2§. 

With the aid of SPEA, we indeed find that the cages of energetically favorable W@Sii2 
and W@Sii4 clusters are simple 3-polytopes with iV= (6,0,2) [Fig. 1(d)] and (3,6,0) [Fig. 1(e)], 
respectively. Varying dcut from I.OA down to 0.25A, we find that ^3 s of both are constants 
(6 and 4; see Table I for the summary of the result). On the other hand, ^3 s of the other 
W@Sii2 and W@Sii4 clusters decrease to 0—2 as dent becomes smaller, whose total energies 
are all higher than those of the (6,0,2) and (3,6,0) cages, respectively. These results support 
the topological argument made above. 

It should be noted that the WQ)Si„ (n~10) clusters may be viewed as metastable phases 
of Si- rich tungsten silicides, which do not appear in equilibrium W-Si phase diagram PB| . 
It is suggested in the previous experiment [^ that the clusters are obtained via a sequential 
growth where a W@Si„ cluster is obtained by attaching additional m Si atoms to a smaller 
W@Si„_m cluster. In order to understand the essence of the growth process on the basis 
of our energetics, we calculate differential binding energy, defined as AE{n)=Eyji^^{n + 
2) - 2EunM + Eundin - 2), where E^Un)= i?t"otT(W) + nE»„tT(Si) - ^^^^(WSiO. A 
negative value of A£'(n) means that generation of a W@Si„ cluster is favorable. We find 
that AE(10)=1.6 eV, AE(12)=-2.8 eV and AE(14)=-0.3 eV. For n=10, stable caging 
of a W atom in simple 3-polytopes [Fig. 1(b)] is unhkely to occur. The cluster should be 
chemically very reactive, because there are too few Si atoms around W to fulfil its reaction 



points. If additional Si atoms can arrive at the cluster prior to its collapse, then subsequent 
growths to W@Sii2 and possibly to W@Sii4 should occur with high probabilities. 

We also find that there exists a threshold cage size for W@Si„ beyond which the cages 
become unstable. In Fig.2, some of calculated structures of W@Sii6 clusters are shown. The 
initial configurations of the cages are taken like simple 3-polytopes with 16 vertices. Upon 
relaxation, all cages we considered are completely distorted and become similar to either 
the (6,0,2) cage plus 2 Si atoms or the (3,6,0) cage plus 4 Si atoms. It appears that n=14 is 
the threshold size in the case of W doping. For a cage beyond a critical size, there are too 
many Si atoms around W to cover its reaction points. 

In conclusion, we show that it is possible to construct Si clusters in fullerene-like 
cages with transition-metal atom doping. A topological discussion suggests that simple 
3-polytopes with maximum numbers of inner diagonals close to the dopant may be good 
candidates of fullerene-like Si cages. First-principles calculation shows that, in the case of 
W doping, W@Sii2 and W@Sii4 are energetically most favorable and also have the cages 
predicted by the above topological picture. 

This work is partly supported by NEDO. 
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TABLES 

TABLE L Topological properties and energetics of W@Si„ clusters. For definitions of n, N, 
^3) ^3 s-iid dcut, see text, di and d2 are the Si-W and Si-Si distances in A, averaged over the 
respective pairs whose separations are less than 2.7A. AE is the total energy ( eV) relative to that 
of the lowest-energy cluster, A^=(6,0,2) for n=12 and A^=(3,6,0) for n=14. 
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14 (3,6,0) 34 4 4 4 4 2.69 2.34 Fig.l (e) 0.00 

14 (4,4,1) 33 11 5 1 2.66 2.35 Fig.l (f) 0.95 

14 (5,2,2) 32 13 4 2 2 2.64 2.39 Fig.l (g) 0.63 

14 (6,0,3) 31 9 7 4 2.63 2.39 Fig.l (h) 0.89 



Figure captions 

Figure 1 (color) Calculated structures of W@Si„ in "fuUerene-like" cages. Double line segments 
are drawn when the inter- atom distance is less than 2.7A. The double line segments 
corresponding to the edges of simple 3-polytopes are colored with blue. In panels (a), 
(g) and (h) , red single hues should be supplemented and also green double lines should 
be removed to retrieve connectivity in the corresponding simple 3-polytopes. See also 
a note |T^ . 



Figure 2 (color) Calculated structures of W@Sii6. Si atoms colored with dark blue are regarded 
"attached" to either the (6,0,2)-cage W@Sii2 (panel (e)) or (3,6,0)-cage WOSiu clus- 
ters (panels (a), (b), (c) and (d)). Other conventions are the same as in Fig.l. 



